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Abstract 

The  well-known  surface  densities  of  bound  charge  o — P • n and  current  K — M x n,  which 
give  rise  to  discontinuities  in  fields  at  a surface,  are  inconsistent  in  their  multipole  order. 

The  electric  quadrupole  contribution  is  missing  from  each.  A consequence  of  including  this 
is  the  appearance  of  a surface  density  of  electric  dipole  moment  P.  Its  effect  on  the  boundary 
conditions  is  derived,  and  a reflection  experiment  on  a bianisotropic  crystal  is  proposed,  which 
in  conjunction  with  a transmission  experiment,  allows  the  V contribution  to  be  measured. 

1.  Introduction 

The  contribution  to  a physical  effect  of  electric  quadrupoles  induced  in  matter  by  long  wavelength 
radiation  is  known  to  compare  in  magnitude  with  that  of  induced  magnetic  dipoles  [1,  2,  3]. 
Consequently,  there  is  an  inconsistency  in  the  surface  densities  of  bound  charge  and  current  [4] 

a = P • n,  K = M x n,  (1) 

where  n is  the  outward  unit  normal  at  a macroscopic  point  on  the  surface  and  P and  M are 
the  electric  and  magnetic  dipole  moments  per  unit  macroscopic  volume.  An  electric  quadrupole 
term  should  be  included  in  K alongside  M,  and  also  one  in  o. 

This  paper  derives  these  two  contributions  and  then  shows  their  effect  on  the  boundary 
conditions  on  the  fields  at  a vacuum-dielectric  interface.  In  particular,  the  electric  dipole  surface 
density  also  emerges  and  its  role  is  examined,  as  well  as  its  contribution  to  a physical  effect. 

2.  Theory  of  Bound  Source  Contributions 

The  retarded  vector  potential  outside  a bounded  distribution  of  currents  in  vacuum  is  [4] 

A{K,t)  = (no/4ir)  f J(r,t-  |(7 1 - r|/c)/|(72.  - r| dv.  (2) 

Jv 

By  expanding  J for  TZ  r,  applying  the  condition  A > U and  averaging  over  a macroscopic 
volume  element  dV,  we  obtain  the  macroscopic  vector  potential  due  to  a bounded  dielectric  of 
volume  V and  surface  area  S.  To  order  electric  quadrupole  (E2)-magnetic  dipole  (Ml)  it  is  [5] 

Aa(R,t)  = (ho/4tt){J  {l/R){Pa  - \Vf}Qap  + ea/?7V^M7)(fF  + J {l/R)(±Qa/3-eai3yMy)npda}. 

(3) 

In  this  R = —71  is  the  displacement  of  dV  from  the  field  point  and  Qap  — X qrarp/dV  is  the 
electric  quadrupole  moment  per  unit  volume.  The  multipole  moment  densities  Pa , Ma,  Qap  are 
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relative  to  an  arbitrary  origin  in  dV  at  the  time  t — R/c.  Similarly  for  the  scalar  potential  [5] 
m ,t)  = l/(47re0){^(l/i?)VQ(-PQ  + lV/3QQ(3)rfF 

+ J [(l/R){(Pa  - ^VpQap)na  + V||P„}  + (P±/J?3)Pj.]da} , (4) 

where  Vy  denotes  differentiation  parallel  to  the  surface  element  da  and 

P a = ~^QapnP  (5) 

is  a contribution  of  E2  order.  From  the  volume  integrals  in  (3)  and  (4),  with  their  l/R  depen- 
dence, the  bound  source  volume  densities  of  current  and  charge  are 

Ja  = Pa  ~ pQap  ~b  capjX  pM^,  p = VQPa  + | XaXpQap.  (6) 

These  satisfy  the  equation  of  continuity  for  bound  sources  V • J = — p.  From  the  l/R  terms  in 
the  surface  integrals  in  (3)  and  (4)  the  surface  densities  of  bound  current  and  charge  are 

Ka  = {\QaP  ~ ea^<yA£y)Tly3,  a = ( P a — ±VpQap)na  + V||7',||.  (7) 

There  still  remains  in  the  surface  integral  in  (4)  the  term  with  a R /R3  dependence.  As  this 
is  characteristic  of  the  potential  due  to  an  electric  dipole,  we  interpret  P in  (4)  as  the  surface 
density  of  electric  dipole  moment  due  to  bound  charge  on  the  surface. 

If  the  bound  sources  in  (6)  are  used  in  the  two  inhomogeneous  Maxwell  equations 


V • E - (pc  + p)/eo,  V x B — + Jc -+- J),  (8) 

where  pc  and  Jc  are  the  corresponding  free  source  terms,  one  obtains  by  comparison  with 

V ■ D = pc,  V x H = D + Jc,  (9) 

the  multipole  expressions  for  D and  H to  order  E2-M1 

Da  — ^Qpa  "b  Pa  f)Qaf)->  Pa  ~ Pa/ 11 0 Ma.  (10) 

The  densities  Pa , Ma,  Qap  are  induced  by  the  wave  fields  and  their  space  and  time  derivatives 
Pai  Pat  VppQ .,  V pEa,  Ba,  Ba , X pBa,  X pBa,  ....  (H) 

Then  to  E2-M1  order 

Pa  = oiapEp  + u~loi'apEp  -F  5aap7V7Ep  -I-  \u;-1a,ap1'V1E0  + GapBp  + u}~lG'apBp,(  12) 
Qa{3  — a^ap  Ery  -f-  Ld  _1 a'-yapEp , (13) 

Ma  = GpaEp—u>  1G'paEp,  (14) 


where  u is  the  angular  frequency.  From  their  definitions  E,  V,  P,  Qap  are  time-even  and  B, 
M time-odd.  Thus  a'Qp,  a'apy , Gap  are  time-odd  and  belong  only  to  magnetic  crystals.  With 
(13)-(14)  in  (10),  constitutive  relations  for  D and  H are  obtained.  However,  these  do  not  satisfy 
Post’s  covariance  requirement  for  a plane  time-harmonic  wave  for  negligible  absorption  [6] 

Da  = AapEp+TapBp,  Ha  = UapEp+XapBp,  Aap  = A*pa,  XQp  = Xpa,  Uap  = -Tpa.  (15) 
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From  V x E = -B  and  with  (10),  (12)-(14)  substituted  into  (9),  one  can  show  that 

Pa  — ocapEp  + u~1[a'Qp\Ep  + $u~l[a'a0y  + a'^a  4-  a^a/3]Vr,Ep 

+ [Gap  ~ k&apG^  — luepjsa'ySalBp  + w 1(G'ap  — kuefasa-f, 5a)Bp,  (16) 

QaP  — ~ku)  iaapy  aPya  "F  a"/ap\Py  = ~~^aPyByi  (17) 

Ma  = [Gpa  — $6apG*y~f  — ^U}CarySa'ySp]Ep  ~ W ( Gpa  — \ti>€a rySay6p)Pp-  (13) 

Terms  in  brackets  [ ] are  time-odd.  From  (10),  (16)— (18)  and  with  E = Eoexp{— i{ujt  — k • r)} 
the  covariant  forms  in  (15)  are  satisfied.  Then  from  (6)  the  Ampere-Maxwell  equation  in  (8)  is 

^aP'y^ pP 7 "h  ^^PO^apBp  — 6 (1®) 

for  a dielectric,  where  eap  is  the  dynamic  permittivity  tensor,  which  to  order  E2-M1  is 

Cap  = eofiaP  + Otap—i&ap+U^kyiAapy—iA'apy),  (20) 

Aap-y  — -tp-iSGaS  ~ ^a^pGps  + \^{aaP^  F aPaj)t  (21) 

Aapy  = —CpySGaS  + CarfSGps  — %U)(aap-y  — apay).  (22) 


We  now  show  that  the  surface  discontinuities  in  (7)  alter  the  usual  Maxwell  boundary  conditions. 


3.  The  Boundary  Conditions 


To  derive  these,  unit  step  functions  are  used  instead  of  the  integral  forms  of  Maxwell’s  equation. 
With  the  surface  element  in  the  xy  plane  and  +z  axis  into  the  medium,  these  functions  and 
their  derivatives  are,  using  the  Dirac  ^-function  and  its  derivative  S' , 


u(z)  = 


( 


1 for  z > 0 
0 for  z < 0 


u{-z) 


( 0 for  z > 0 
( 1 for  z < 0 


du{±z ) 
dz 


= ±«(*), 


£^(±2) 

dz 2 


= ±5'(*),  (23) 


The  total  bound  current  and  charge  densities  are  then 


Jt(R)  = ti(*)Ji(R)  +u{-z)  J2(R)  + tf(z)K(r),  (24) 

Pt(R)  = u(z)pi(R)  + u(-z)p2(R)  4-  S{z)(7{t)  + S'(z)z  ■ "P(r),  (25) 

where  medium  2 is  the  vacuum,  n = — z,  and  r lies  in  the  xy  plane.  Using  (23)-(25)  in  V- J = —p 
and  equating  the  coefficients  of  S(z)  one  obtains  cr  in  (7)  and  of  S'(z)  one  finds 

Kz  = -Vz  = \QZZ,  (26) 


which  is  K in  (7)  for  a = z.  Since  these  two  results  confirm  (24)  and  (25),  we  similarly  take 


E(R)  = u(z)E1(R)  + ti(-z)E2(R)  + 6(z)£  (r),  (27) 

B(R)  = ti(«)Bi(R)  + «(-*)B2(R)  + 5(z)e(r),  (28) 


where  £ (r)  and  S(r)  are  surface  fields.  With  these  in  (8)  one  obtains  the  boundary  conditions 

E\x  ~ Bzx  — Vx5z,  Ely  — E%y  = X?y£z,  E\z  — E2z  — ((J  — tlaz^Ea^z)  Aoi  (29) 

B\x  ~ B2x  = yo (Ay  4-  ti yZ£z),  ( B\y  — B%y)  — — po (Kx  + t\xz£z)i  B\z  = B2Z , (30) 

&z  = Pzl^lzz  = -\Qzzhlzz-  (31) 

Thus  to  order  E2-M1  the  tangential  components  of  E are  no  longer  continuous  across  a boundary, 
whereas  the  normal  component  of  B is. 
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4.  The  Role  of  the  Surface  Density  of  Electric  Dipole  Moment 

As  is  evident  from  (29)-(31)  the  boundary  conditions  reduce  to  the  Maxwell  forms  when  "P  (or 
Qap)  is  zero.  In  the  electric  dipole  (El)  order  Qap  = 0.  Thus  the  Maxwell  conditions  are  an 
approximation,  applying  only  in  the  El  order.  Prom  (5)  and  (17),  "P  exists  only  for  magnetic 
media.  Also,  the  tensor  SQja7  in  (17)  may  vanish  for  symmetry  reasons,  as  for  propagation  along 
the  main  axis  of  antiferromagnetic  Cr203,  even  though  it  possesses  the  time-odd  tensor  a'a^. 

5.  Application 

The  simplest  magnetic  crystal  is  an  antiferromagnet  (since  a'ap  — 0)  that  is  centrosymmetric 
(since  its  time-even  odd-rank  polar  tensor  aa ^7  and  even- rank  axial  tensor  G'ap  vanish).  The  only 
such  symmetry  for  which  an  effect  exists  at  normal  incidence  is  4/mmm.  Its  Fresnel  reflection 
matrix  for  normal  incidence  parallel  to  its  C2  axis  can  be  shown  to  be  [7] 

/?  = f ~ l)/(no  + *)  rP* 

[ rsp  ~(ne  - 1 )/(ne  + 1) 

where  n0  and  ne  are  the  ordinary  and  extraordinary  refractive  indices, 

rsp  = rps  = noc[Kne/ (n0  + ne)  - Si23]/(n0  + l)(ne  + 1),  K = 2 Gn  + w(a'123  + a312),  (33) 

and  S'i23  is  the  only  component  that  exists  of  the  tensor  in  (17)  which,  because  of  (5),  is  the 
surface  electric  dipole  term.  The  components  in  (33)  are  relative  to  crystallographic  axes.  The 
matrix  in  (32)  is  identical  in  form  to  that  of  (>203  for  the  same  geometry.  The  effect  in  Cr203  is 
a rotation  of  the  plane  of  linearly  polarized  light,  which  has  previously  been  measured,  yielding 
an  experimental  value  of  rsp.  Thus  rsp  should  be  measurable  for  a 4/mmm  crystal.  A different 
combination  of  K and  S123  enters  the  birefringence  in  transmission 


(32) 


n+  - n_  = c/i0(35i23  - K).  (34) 

Prom  these  reflection  and  transmission  experiments  a value  can  be  obtained  of  6123  and  hence  of 
the  contribution  of  the  surface  density  of  electric  dipole  moment  "P  to  the  boundary  conditions. 
This  would  provide  a test  of  our  theory. 
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